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ABSTRACT 

We present an interferometric study of the continuum surface of the red supergiant star Betelgeuse 
at 11.15 /.tm wavelength, using data obtained with the Berkeley Infrared Spatial Interferometer each 
year between 2006 and 2010. These data allow an investigation of an optically thick layer within 
1.4 stellar radii of the photosphere. The layer has an optical depth of ~1 at 11.15/im, and varies 
in temperature between 1900 K and 2800 K and in outer radius between 1.16 and 1.36 stellar radii. 
Electron-hydrogen atom collisions contribute significantly to the opacity of the layer. The layer has 
a non-uniform intensity distribution that changes between observing epochs. These results indicate 
that large-scale surface convective activity strongly influences the dynamics of the inner atmosphere 
of Betelgeuse, and mass-loss processes. 

Subject headings: infrared: stars — stars: individual (Betelgeuse) — circumstellar matter — stars: 
mass- loss — techniques: interferometric 



1. INTRODUCTION 

Mass-loss from massive stars in the red supergiant 
phase is fundament al to th e chemical enrichment of the 
interstellar medium (jGehrdllQSS ). Red supergiants, with 
initial masses of 10— 50M (T). are known to lose mass at 
rates of up to lO'^^Mrr^/vr (|Humphrevill986D . Mass-loss 
rates are commonly estimated indirectly, using observa- 
tions of infrared continuum excesses fe.g.. ^S tencel et all 
Il989f ) or measure ments of C O rotational line intensi- 
ties (e.g.. Ide Beck et aD 120101 ). Outflowing material is 
also directly surveyed through observations of system- 
atic s pectral line displacements (jAdams fc MacCormackl 
I1935D and interfero metric studies over multiple epochs 
()Bester et aI1ll996( ). The physical processes involved in 
red supergiant mass-loss are however not well under- 
stood. While mass-loss from less massive stars, such 
as AGB stars, is understood in terms of significant ra- 
dial pulsations and winds driven by radiation pressure 
on dust grains, such mechanism s are not applicable to 
red supergiants (llosselin & Plez 2007") , except those in 
a superwind phase (Yoon & Cantiello 2010). Rotation 
in stars that are spun-up as they condense from pro- 
tostellar clouds could play a role in driving mass-loss 
(jHeger fc Langed 119981). as could Alfven and acoustic 
waves (e.g.. iHartmann fc MacGregoilll980[ ). In the case 
of Betelgeuse (a Orionis, spectral type M2Iab), th e el- 
evation of cool photospheric gas (Lim et al.' ' 1998D and 
gas a t chromospheric temperatures (Gilliland fc Dupred 
Il996f l has be en attributed to large-scale surface con vec- 
tive activity ()Schwarzschildlll975l : lAntia et al.lll984l ). 
Betelgeuse is perhaps the archety pal red supergiant. 



Its la rge angular diameter (0.047", iMichelson fc Peaii 
|1921| ) and high apparent magnitude {Mk = —4.38) make 
it an excellent subject for high spatial resolution stud- 
ies. The continuum angular diameter of Betelgeuse has 
been measured u sing interferometri c techniques at radio 
(|Lim et al.lll998D . mid-IR f8-12 ^m. lWeiner et al.ll2003al : 
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and ultraviolet (|Gilliland fc Dupreelll996[) wavelengths, 
and significant variation of the diameter with wave- 
length is observed. The apparent mid-IR diameter is 
also observed to vary significantly on timescales of years 
(jTownes et al.l 120091 ). Numerous interferometric cam- 
paigns, from optical to mid-IR wavelengths, have further 
revealed up to three hotspot s on th e stellar surfac e (e.g., 
[Busc her et al.. 1990 : Youn g et all UDOO; Haubois eFaLl 
120091: iTatebe et all 120071 ). 'These spots range in inten- 
sity from a few to a few tens of percent of the intensity 
of the stellar disk. The detection of these spots is gen- 
erally interpreted as direct evidence for the presence of 
a few giant convectio n cells on the stellar surface (e.g., 
iChiavassa et"aL|[20To[ ). 

A profusion of circumstellar components have been 
identified surrounding Betelgeuse. Furthest from the 
star, an analysis of mid-IR. spectra and interferometric 
data bv lDanchi et al.l (|1994f ) showed the presence of two 
thin shells of dust at approximately 40 and 80 stellar 
radii (i?*), with temperatures of ^400 K and ^200 K re- 
spectively. Later work bv lBester et al.l (|1996[ ) revealed a 
third, more compact dust shell at 4i?* that was not ap- 
parent in previous data. The dust condensation zone and 
acceleration mechan ism are not ident ified. Radio contin- 
uum observations bv lLim et al.l (|1998f ) at wavelengths be- 
tween 0.7 cm and 6 cm, probing spatial scales of 2 — 7i?* , 
revealed the presence of neutral gas at temperatures be- 
tween 3500 K and 1300 K. These temperatures are con- 
sistent with the temperature profile of the Betelgeuse 
atmosphere in these regions, as derived from mean tem - 
perature and density models (|Harper et al.|[200ll l2009f) . 
At 2i?*, the gas distribution is asymmetric, indicating 
an anisotropic mechanism fo r elevating the gas above 
the phot osphere. Ultraviolet (jCilliland fc Dupreill996() 
and Ha (jHebden et al.lll987D imaging observations have 
shown that this cool gas co-exists with asymmetrically- 
distributed hot chromospheric material. Observations of 
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[Fe II] lines by iHarper et al.l (|2009[ ) that correspond to 
the dominant cool atmospheric component showed, sur- 
prisingly, that this material was at rest with respect to 
the star. The circumstellar envelope as viewed across 
the JHK wavelength bands was imaged using a 'lucky- 
imaging' approach c ombined with adaptive optics by 
iKervella et all (|2009[) . and was found to include a plume 
extending to 6i?* . 

Closer to the star, strong evidence exists for the 
presence of a cool molecular layer or shell directly 
above the photosp here within 1.5^* (e.g. lOhnaka 200l; 
iPerrin et al.l2007j ) . This layer is optically thick at mid-IR 
wavelengths but is optically thin in the near-IR. Spectral 
signatures of water from Betelgeuse and the red super- 
giant /i Cephei, found using data from th e Infr ared Space 
Observatory (ISO), were interpreted bv iTsuj i (2000a, b) 
as originating from a 1500 K layer surrounding these 
sta rs. Modeling of nea r- and mid-IR interferometric data 
by IPerrin et al.l ()2004|) confirmed the presence of such a 
layer, with a mid-IR optical depth of ~2 and a temper- 
ature of ~2000K. Further work bv lTsuiil (|2006[) combin- 
ing both interferometric and spectroscopic data placed 
the layer surrounding Betelgeuse at beyond 1.3i?*, with 
a temperature of 2250 K. T he layer was shown not to be 
solely composed of water bv lVerhoelst et all ()2006D . who 
suggested the presence of amorphous alumina (AI2O3) to 
explain the high mid-IR opacity. The spectrum of wa- 
ter lines from Betelgeuse is indeed more complex than 
revealed by the ISO data; signatures o f water at photo- 
spheri c tem peratures we r e dete cted bv lJennings &: Sadal 
(|1998D and iRvde et al.l (|2006D . Most recently, wide- 
band, spect rally dispersed mid- IR interferometric obser- 
vations by IPerrin et al.l ()2007f) were fitted by a layer 
containing SiO, amorphous alumina and water. High 
resoluti on measurem ents of the CO overtone band near 
2.3 um ([Ohnaka et al. 2009) and one-dimensional imag- 
ing (jOhnaka et al .11 2 Oil.) were modeled with a large patch 
of gas near 1.45i?* with motions of 10— 15kms~^ with re- 
spect to the star. This emerging model of a photosphere 
and overlying layer for a selection of red supergiant stars 
allows a better understanding of dust formation and re d 
supergiant wind acceleration fe.g. lVerhoelst et aI]|2009D . 
leading towards a convergent picture of red supergiant 
mass loss. 

Studies of the extended Betelgeuse atmosphere form a 
basis for a better understanding of the atmospheric dy- 
namics of red supergiant stars. The existence of a cool 
layer with high opacity at mid-IR wavelengths directly 
above the photosphere is now established for Betelgeuse. 
The size, temperature and composition of the layer have 
been investigated, assuming spherical symmetry. Little, 
however, is known about the spatial structure and time- 
variability of this layer. Consistent instrumentation and 
modeling procedures have also not been applied over mul- 
tiple observation epochs. In this paper, we present obser- 
vations of the mid-IR continuum of Betelgeuse usi ng the 
Berkeley Infrared Spatial Interferometer (ISI, Hal e et aD 
12000^. Interferometric data using three telescopes were 
gathered during each year between 2006 and 2010, with 
angular resolutions sufficient to model the stellar shape 
and potential large-scale features. We describe the ob- 
servations and image-modeling processes in §2, and per- 
form fits to a pliotosphere-and-layer model for our ob- 



TABLE 1 
Observing log. 



Year 




Dates (UT) 


11.15 flux density (10'^ Jy) 


2006 




8, 9, 10 Nov; 7 Dec 


4.2±0.2 


2007 




14, 15, 16 Nov 


4.1±0.3 


2008 




22, 23, 24 Sep 


3.8±0.5 


2009 


4, 


5, 7, 9, 18, 20, 22 Nov 


4.4±0.1 


2010 


5, 


17, 18 Nov; 9, 13 Dec 


4.9±0.3 



servations in §3. In §4, we discuss the nature of the 
observed non-uniformities in our model images. We ex- 
amine the implications of our results for models of the 
Betelgeuse atmosphere and red supergiant mass-loss in 
§5, and present our conclusions in §6. 

2. ISI OBSERVATIONS OF BETELGEUSE 

2.1. Observational methods 

The ISI consists of three 1.65 m aperture telescopes 
with a heterodyne detection system at each telescope. 
The double-sideband detection systems are sensitive to 
9 — 12 /xm radiation within single sideband widths of 
~2.7GHz. The ^'^C02 laser local oscillators can operate 
at numerous wavelengths corresponding to vibrational 
molecular transitions between wavelengths of 9 /im and 
12 /im. For the observations reported here, the local os- 
cillators at each telescope were tuned to a wavelength of 
11.15 /.tm, chosen so as to exclude strong s tellar and tel- 
luric sp ectral lines from the detection band (jWeiner et al.l 



The van Cittert-Zernicke theorem equates the com- 
plex cross-correlation of the electric fields detected at two 
spatially-separated telescopes with the two-dimensional 
Fourier transform of the source brightness distribution 
within the telescope fields of view. Instead of measuring 
complex cross-correlations, the ISI, like all optical/IR in- 
terferometers, measures the normalized cross-correlation 
amplitude, known as the Michelson fringe visibility. The 
Michelson fringe visibility, V , is given by 



V = 



P. 



fringe 



(1) 



where the fringe power, Pfrmge is the cross-correlation 
amplitude between signals measured at two telescopes, 
A and i?, and Ia and Ib are the source powers mea- 
sured at telescopes A and B respectively. Visibility am- 
plitudes are recorded for each ISI telescope baseline pair. 
The visibility phase in the mid-IR for a given baseline 
is, however, strongly affected by randomly varying at- 
mospheric pathlengths above each telescope, and cannot 
be recovered from the data. The ISI instead measures 
the closure phase, or the sum of the measured visibility 
phases for each baseline. The closure phase is indepen- 
dent of opt ical pathlength differences between pairs of 
telescopes (jJennisonl I1958D and it provides information 
about the point-inversion symmetry of the source. The 
spatial frequency sampled by the projection of a baseline 
of length 6, perpendicular to the direction of a source, 
at a wavelength A, is proportional to b/X. By obtaining 
measurements at different source hour-angles, multiple 
spatial frequencies at different position angles are sam- 
pled. 

Visibility and closure phase data were recorded for 
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Fig. 1. — Plots of the ISI measurements of Betelgeuse from 2006—2010. The left panels show visibilities, the center panels show closure 
phases, and the right panels show UV coverages, for each observing epoch (2006 to 2010 from top to bottom). The visibilities for each 
baseline, denoted 1-2, 2-3 and 3-1, are depicted by diamonds, triangles and squares respectively (note that the spatial frequency ranges 
vary from year to year). The fitted model visibilities (Table 2) are shown as red, green and magenta lines for baselines 1-2, 2-3 and 3-1 
respectively, and as magenta lines for the closure phases. The closure phases are plotted against the position angle. East of North, for 
baseline 2-3. In the plots of the UV coverage, baseline 1-2 is represented by red diamonds, baseline 2-3 is represented by green triangles, 
and baseline 3-1 is represented by blue squares. 
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Betelgeuse over timespans of five weeks or less each year 
between 2006 and 20f0 (see Table f). During this time 
period, the three ISI telescopes were configured in an ap- 
proximately equilateral triangle with baseline lengths be- 
tween 34 m and 40 m. For Betelgeuse, this configuration 
allows spatial frequencies between 20 and 37SFU (Spa- 
tial Frequency Unit; f SFU= f x fO^ cycles radian"^) to 
be sampled. For all observations, tip-tilt corrections were 
applied to center the source on .ftT-band guider cameras 
operating at frame rates of 27 Hz. A standard observ- 
ing sequence consisted of 4.5-minute dataset intervals of 
staring at the source and position-switching between the 
source and blank fields on either side of the source. The 
position-switching measurements moved the telescopes 
between the source and the sky every 15 s. The source 
powers were measured using additional 140 Hz chopping 
between the source, or the sky, and cold loads with tem- 
peratures near that of the sky. A single visibility mea- 
surement consisted of either: (1) a fringe power measured 
during a chopping/position-switching dataset, normal- 
ized using the source powers during that dataset; or (2) 
a fringe power measured during a staring dataset, nor- 
malized by weighted averages of source powers measured 
during the nearest chopping/position-switching datasets. 
Closure phases were calculated for each triplet of visibil- 
ity measurements. 

Visibility magnitudes and closure phase zero-point off- 
sets were calibrated using observations of Aldebaran (a 
Taurus) conducted on the same night. Aldebaran was 
assumed to have a diameter of 20 mas and zero closure 
phase. The visibility calibrations for the two kinds of 
visibility measurements were different because of the dif- 
ferent times spent on-source during chopping and staring 
datasets. Data that were alfected by systematic errors or 
poor atmospheric conditions were discarded. The total 
flux densities of Betelgeuse during each observing epoch 
were also estimated using Aldebar an, assumed to have a 
flux density at 11.15 ^m of 615 Jy (|Monnier et al.lll998f ). 
as a calibrator. These total flux density values are given 
in Table 1, and represent the total power of Betelgeuse 
within the 5"x5" fields of view of the ISI detectors. 



2.2. Image modeling 

Image models were fitted to the calibrated visibil- 
ity and closure phase datasets for each epoch. The 
data, analytic fits, and coverage of the UV spatial fre- 
quency plane, are shown in Figure 1. Effective model- 
independent image reconstruction was not possible be- 
cause of the sparsity of the UV coverage in all epochs. 
Our model images include a centered uniform disk (UD), 
and up to two point sources offset from the image cen- 
ters. The uniform disks represent the apparent stellar 
surface, and the point sources represent localised inten- 
sity fluctuations in this surface. Emission from the dust 
surrounding the star at large angular scales, while con- 
tributing to the total detected intensity, was assumed 
not to contribute to the visibility at the sampled spatial 
frequencies. The intensity of each image component was 
fitted to a fraction of the total intensity. Other free pa- 
rameters included the UD radius and the positions of the 
point sources. The model complex visibility at a point 



(u, v) in the UV plane was given by 



, 2AJi(27rrVu2 + t,2) 

V[U, V) = r^=^= h 

27rrv + 

p^^-2iTi{uxi+vyi) _|_ p^^-2iTi(uX2+vy2) 

where A is the fraction of the total intensity contributed 
by the UD, r is the uniform disk angular radius, Ji de- 
notes a Bessel function of the first kind and of order 
unity, P„ represents the fractions of the total intensity 
of the point sources, and a;„ and ?;„ represent angular 
offsets to the West and North for the point sources, with 
n = {l,2}. 

We used a weighted least-squares fitting technique for 
each epoch, with each baseline and the closure phase data 
given equal weights. The differing quality of data from 
different epochs necessitated using models with different 
numbers of free parameters to fit different epochs. Data 
from 2006 were fitted with a single point source in addi- 
tion to a UD, data from 2007 and 2008 were only fitted 
with UDs, and data from 2009 and 2010 were each fitted 
with a UD and two point sources. An additional point 
source was added to the fit for a given epoch only if the 
reduced values calculated for the full dataset, and for 
the closure phase data alone, were both decreased upon 
the addition of the point source. We also attempted to 
fit other image models. These included uniform ellipses, 
uniform ellipses with point sources, and UDs with offset 
uniform disks or offset two-dimensional Gaussian profiles 
in place of the point sources. None of these trials, how- 
ever, converged to a satisfactory image, and the angular 
sizes of the fitted point sources are not determined. The 
fractions of the total intensity contributed by the point 
sources were not constrained to be positive, as our data 
are sensitive to both bright and dark features on the stel- 
lar disk. 

Simulations were conducted to assess possible degen- 
eracies in fitted image model parameters given the spar- 
sity of the UV coverage of our datasets, in order to ensure 
that our data were not being over-fitted. No significant 
degeneracies were found for any epoch. For example, 
the 2006, 2009, and 2010 datasets constrain both the lo- 
cations and brightnesses of the point sources. Closure 
phase measurements at multiple array position angles 
constrain the position angles of the point sources, and 
both the closure phase measurements and the visibility 
measurements constrain the angular offsets of the point 
sources from the image centers. 

The results of the fits, as well as the reduced for each 
epoch, are given in Table 2. All the fits have reduced x^ 
values that are greater than unity. This implies that the 
data contain minor features besides Gaussian noise that 
are not being fitted by the free parameters of the image 
models. 

ISI data fo r Betelgeuse obtaine d during 2006 were also 
analyzed by iTatebe et al.l (|2007| ) . Despite our choice of 
different observing dates, a different method of calibra- 
tion and slightly differe nt image fitti n g pro cedures, our 
results match those of iTatebe et al.l ()2007[ ) within the 
margins of error. We however do not adopt Tatebe et 
al.'s results for the 2006 epoch in order to maintain the 
same analysis procedure across all epochs. 
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TABLE 2 

UD AND POINT SOURCE(s) MODELS FITTED TO VISIBILITY AND CLOSURE PHASE DATA. 



Year UD fraction UD radius Point source fraction Pt src x (mas) Pt src y (mas) Reduced 



2006 
2007 
2008 
2009 

2010 



0.51±0.02 
0.54±0.01 
0.48±0.02 
0.557±0.008 

0.61±0.03 



24.5±0.8 
24.8±0.3 
24.7±0.8 
26.4±0.2 

26.4±0.3 



0.04±0.01 



0.009±0.001 
0.007±0.001 
0.025±0.009 
0.009±0.004 



-2.4±0.1 



-25.0±0.3 
18.3±0.3 

-29.7±0.3 
8.3±0.3 



-23.7±0.2 



10.5±0.2 
-20.2±0.2 
-9.3±0.1 
-25.9±0.1 



2.683 
2.053 
2.416 
4.565 

3.65 



2.3. The changing angular diameter of Betelgeuse 

Using ISI observa tions over the period 1993—2009, 
iTownes et al.l ()2009D reported a systematic decrease of 
15% in the mid-IR continuum UD diameter of Betel- 
geuse. Here, we augment this finding with new mea- 
surements of the UD diameter of Betelg euse. Figure 
2 show s the data plotted in Figure 1 of iTownes et al.l 
(p009t ). with the new results also included. Though the 
present wo rk re-analyzes some previous observations ex- 
amined by iTatebe et al.l ()20Q7f ) , the diameter measure- 
ments are obtained fro m different data sets a nd analyses 
from those reported bv To wnes et al.l (|2009D . The mea- 
sur ement of the UD di ameter of Betelgeuse at 11.15 fj,m 
bv IPerrin et al.l (|2007[ ). obtained using the Very Large 
Telesc ope Interferorneter in 2003, is also plotted. Re- 
cently, lOhnaka et al.l ()201lD showed that K-band diame- 
ter measurements, made over a similar time period, vary 
less than the mid-IR diameters plotted in Figure 2. 

While our results from 2006—2008 are consistent to 
within the ma rgins of error with the trend towards 
a smaller size ()Townes et al.l I2OO90 , the UD diameters 
of Betelgeuse in 2009 and 2010 are significantly larger 
than those during 2006—2008. Our analysis method for 
the 2007 and 2008 datasets essentially matched that of 
ITownes et al.l (f2009i ) , in that only UDs were fitted to the 
data. Our technique however differed somewhat in the 
methods of calibration and fitting. We consider the vari- 
ability in the apparent size of Betelgeuse in the mid-IR 
further in §5.2. 

3. A PHOTOSPHERE-AND-LAYER MODEL FOR 
BETELGEUSE 

3.1. The effective temperature 

Effective temperatures of stellar surfaces are among th e 
fundamental parameters of stars (jvan Belle et al.ll20Tol) . 
The ISI measurements of the flux density of Betelgeuse, 
and of the UD sizes and fractions of the total intensity, 
allow the effective temperature of the 11.15 /im apparent 
stellar surface to be determined at each observing epoch. 
We use a Planck function for the power emitted by a 
blackbody to express the effective temperature, T^ff, as 

where h is Planck's constant, v is the observing fre- 
quency, fcs is Boltzmann's constant, fl is the solid an- 
gle subtended by the stellar disk, c is the vacuum speed 
of light, and P is the measured power from the stellar 
disk in units of Wm~^ Hz"^. We do not use the widely- 
applied technique of calculating effective temperatures 
using estimates for the bolometric flux from the star (e.g.. 



IPerrin et alll2004l: IBovaiian et"al]|2009[ ) for two reasons: 
radiation at mid-IR wavelengths is strongly affected by 
thermal emission from circumstellar dust shells and bolo- 
metric correction factors at these wavelengths are unique 
to individual stars; and also because the correspondance 
between the 11.15 /xm scattering surface and the photo- 
sphere of Betelgeuse is not certain. 

The calculated effective temperatures for the observed 
stellar disk of Betelgeuse during each observing epoch are 
plotted in Figure 3. We do not account for attenuation 
by the interstellar medium, estimated as modifying the 
obser ved flux density by a factor of 0.975 (jBester et al.l 
I1996D . We also do not account for absorption of the stel- 
lar radiation by the dust that surrounds Betelgeuse at 
angular scales of 1" or larger, or for the emission from 
this dust in the line of sight to the stellar face. This 
is because the dust shells are extremely optically thin 
at mid-IR wavelengths when viewed radially, with a to- 
tal line-of-s ight optical depth of approximately 0.0065 
(iDanchi et al. 1994). 

Our measured effective temperatures of the mid- 
IR surface of Betelgeuse are consistently lower t han 
the value of 3641 K derived by IPerrin et afl (|2004D for 
the photospheric component of a photosphere-and-layer 
model. Furthermore, significant variability is evident be- 
tween the effective temperatures measured at different 
epochs. We therefore consider sources of mid-IR contin- 
uum opacity in the atmosphere of Betelgeuse in an at- 
tempt to elucidate the scattering surface probed by our 
observations and the cause of the observed variability. 

3.2. The 11.15 fim scattering surface of Betelgeuse 

Stellar atmospheres have opacities that vary with 
wavelength. Stars appear smaller at wavelengths where 
the atmospheric opacity is smaller. The surface of a 
star, when observed at a given wavelength, can be ap- 
proximated as the location where the optical depth at 
that wavelength, measured from the observer, is approx- 
imately unity. The photosphere of a star, however, is 
the apparent stellar surface observed at a wavelength 
where the opacit y matches the Rosseland mean opacity 
(|Rosselandl I1924D . The Rosseland mean opacity, kr, is 
given by 

where B{i') is the intensity of radiation from a point in 
a stellar atmosphere, T is the temperature at that point, 
and k{i>) is the opacity at a frequency v. Regions of a 
stellar spectrum with large radiation intensity but low 
opacity contribute most to determining the Rosseland 
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Fig. 2. — The UD diameter of Betelgeus e at 11.15 um measured using the ISI over the period 1993—2010. Data between 1993 and 2009, 
plotted as diamonds, are taken from Tow nes et ahl 12009^. and additional data between 2006 and 2010, plotted as squares, are from the 
present work. The 11.15 /im UD diameter measurement of iPerrin et al.i (120071 ). obtained in 2003 with the VLTI, is plotted as a triangle. 
The Modified Juhan Day (MJD) corresponding to a given Juhan Day (JD) is MJD=JD-2400000.5. 
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Fig. 3. — The effective temperature of the apparent 11.15 fim stellar disk of Betelgeuse for each ISI observing epoch in the present work. 
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mean opacity. iWeiner et al.l (|2003aD concluded that the 
mid-IR continuum opacity of a gas with solar atomic 
abundances under surface conditions present on red giant 
stars closely matches the Rosseland mean opacity of such 
a gas. 

Molecular layers surrounding red giant stars have, how- 
ever, been shown to significantly change the mid-IR opac- 
ity from that of a solar abundance atomic gas, while re- 
mainin g relatively transparent in s horte r wavelength IR 
bands (Wcincr 2004: 'Perrin^eLaLl 12004: 'Verhoelst ct alj 
Il006; Woodruff ct al. 2009). TheTPcrrin et al. (200j 
model for Betelgeuse, combining K-han d visibility dat a 
and a 11.15 /im ISI visibility curve (.Weiner et al.ll2000l) . 
includes a shell with a -fiT-band optic al depth of . 06 and 
a 11.15 ^m optical depth of 2.33. iPerrin et afl (|2007[ ) 
showed that this shell was composed of water, amorphous 
alumina dust and SiO gas, and that the 11.15 /im opacity 
was dominated by alumina. The presence of alumina in 
a layer s urrounding Betelgeuse , within 1.5i?*, was sug- 
gested by I Verhoelst et all ()2006f ) as the only dust species 
that could condense at the temperatures present in this 
region. Alumina dust grains forming close to stars pos- 
sibly play a role as nucleation sites f or dust species that 
condense at lower temperatures (e.g. lSalpeterill977l ). 

Bremsstrahlung interactions between electrons and 
neutral hydrogen atoms also contrib ute to the mid-IR 
conti nuum opacity in red giant stars ()Tatebe &: TownesI 
200fih. This e missio n mechanism was adduced by 
Reid & MentenI ()1997l ) to explain the variation with 
wavelength of the apparent sizes of the radio photo- 
spheres of long-period varia bles, and ascribed to Betel- 
geuse by iLim et all ()1998[ ). We discuss this opacity 
model further in §3.4. 

3.3. Fits to a photosphere- and-layer model 

Prior work, the dramatic variations in apparent mid- 
IR diameters of Betelgeuse, and the low effective tem- 
peratures all provide clear motivation to model the new 
ISI observational results in terms of a photosphere and 
a molecular layer . Our method is similar to that of 
IPerrin et aP ([2007^ in that we modeled a spherically sym- 
metric layer surrounding an opaque stellar photosphere 
with an assumed diameter of .04371" and an as sumed 
surface temperature of 3641 K (IPerrin et al.l 1200^ . Vis- 
ibility and closure phase models for any observed point 
sources were subtracted from the data during each ob- 
servation epoch in order to allow the model to be spheri- 
cally symmetric. We did not attempt to fit an asymmet- 
ric layer because of the large number of free parameters 
involved. The free parameters of our model were: the 
temperature, T, of the layer; the optical depth, r, of the 
layer along the line of sight to the center of the star; and, 
the outer radius, R, of the layer. 

The layer was assumed to have uniform temperature 
and opacity, and was also assumed to have an inner ra- 
dius that corresponds to the photospheric radius. Radia- 
tion within the ISI band was assumed to arise only from 
blackbody emission by the layer. 

The modeling procedure involved fitting visibilities de- 
rived from a model image, created using a trial set of free 
parameters, to the measured visibilities. The measured 
visibilities for each epoch, and their errors, were scaled 
by the measurements of the measured power emitted by 
the UD surface of Betelgeuse. The flux density of a pixel 



TABLE 3 

Temperature, optical depth and outer 

RADIUS OF a layer FITTED TO VISIBILITY 
AND FLUX DENSITY MEASUREMENTS. 



Year 


T (KxlO^) 


T 


R (mas) 


2006 


1.9±0.7 


0.8±0.1 


28.6±0.7 


2007 


2.5±0.9 


0.9±0.1 


25.5±0.9 


2008 


2±1 


1.2±0.2 


25±1 


2009 


2.3±0.8 


0.99±0.06 


26.8±0.8 


2010 


2.8±0.6 


1.05±0.03 


29.8±0.6 



in the model image grid was calculated as: 
for z < R^, 

I = S{v, T,)f7e-"^^(^) + S{v, T)Vt{l - e-"^-^(^)), (4) 

and, for z > R^,, 

I = S{iy, T)n{l - e-2"si(2)). (5) 

Here, L{z) describes the optical path length through the 
layer, is the absorption coefficient, z is the radial dis- 
tance of a pixel from the image center, n is the solid 
angle subtended by a single pixel, S{u, T) is the specific 
intensity radiated by a blackbody of temperature T at 
frequency v, and = 3641 K is the assumed tempera- 
ture of the photosphere. L{z) and as are given by 

L{z) = Vi?' - z2 _ ^^2 _ 2 < i?, (6) 

i(z) = Vi?^-z2, z>i?, (7) 
as^T/{R~R,). (8) 

The pixel size of the model images was 1.2 x 1.2 milli- 
arcseconds. For each iteration of the fit, a model image 
was generated and integrated along one dimension, and 
a one-dimensional Fourier transform was performed in 
order to generate a set of visibilities. Each model vis- 
ibility was compared to measured visibilities that were 
closest to it (i.e., within ±0.5 SFU) in spatial frequency. 
A Levenberg-Marquardt fitting algorithm was employed 
to calculate the free parameter set that minimized the 
global x^- 

The results of the fit for each observing epoch are given 
in Table 3. Errors for each parameter were calculated us- 
ing the resulting covariance matrix. In order to confirm 
the validity of our analyses, we added the visibility and 
closure phase models of the point sources, that were ear- 
lier subtracted from the visibility data from each epoch, 
to the best-fit model visibilities, values relating this 
sum to the original visibility data were calculated for each 
epoch, and were found to be comparable to the values 
of the fits to UDs and point sources shown in Table 2. 

The changes with time of the fitted layer tempera- 
tures and outer radii are similar to the changes with 
time in the effective stellar temperatures and the appar- 
ent UD radii respectively. All our fitted layer radii and 
temperatures are roughly consistent with the 1.3— 1.5-R* 
(0.057" to 0.065"diameter) and 1500-2000 K layers mod- 
eled by most autho rs (Tsuii 2000a: Ohnaka 2(^1 iTsuiil 
I2OOI iVerhoelst et"al.u2006; .Perrin et al. .2007.) . The val- 
ues we calculate for the 11.15 /itm optical de pth of ^1 are 
half of those resulting from the models of IPerrin et al.l 
p004.) and .Perrin et al. (2007,). 
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3.4. Modeling the layer opacity 

The lack of wideband mid-IR spectral information dur- 
ing our observing epochs limits modeling of molecular 
gas and dust species that contribute to the layer opacity. 
We are able, however, to investigate some basic physi- 
cal characteristics of the layer using t he op acity mecha- 
nism described bv lTatebe fc TownesI (|2006D . This work 
showed that, for solar abundance gas at temperatures 
similar to those in red giant atmospheres, collision pro- 
cesses between neutral hydrogen atoms and free elec- 
trons dominate the mid-IR opacity. For temperatures of 
2500 K to 3000 K, the primary sources of free electrons 
were found to be thermally ionized Na and Al. The abun- 
dances of H~, H2 and H+ relative to H were found to be 
too low to contribute to the opacity through collisions 
with electrons, and photoionization and photodissocia- 
tion processes were also found to be insignificant. The 
absorption coefficient due to electron- hydrogen atom col- 
lisions, first derived bv lDalgarno fc Land p966l) . is given 
in simplified form as 

kBTneUHA 
cte = 5 (9) 

where rig is the electron number density. Tin is the hydro- 
gen atom number densit y, and A is expressed as a third- 
order polynomial in T (jTatebe fc TownesI I200I . This 
opacity mechanism is used to expl ain the sizes of the ra- 
dio photospheres of red giant stars (|Reid fc Mentenlfl997l : 
Lim et al. 1998), at different wavelengths. For example, 
Lim et al.l (|1998l) found that the radio photosphere of 
Betelgeuse appeared larger, yet cooler, for longer wave- 
lengths. 

We consider the implications of this mechanism domi- 
nating the mid-IR continuum opacity. Equation 9 can be 
used to calculate hydrogen atom number densities from 
the fitted layer optical depths. The ele ctron number den- 
sities were obtained from Table 2 of iTatebe fc TownesI 
(|200l) . which gives the fractional ionizations, relative to 
njj, of various metals at different temperatures. We fur- 
ther used the hydrogen number densities to calculate the 
total layer mass, Miayen for each epoch, using a distance 
of 197±45pc to Betelgeuse (jHarper et al.ll20Q8D . Our re- 
sults are given in Table 4. 

We find masses for the layer of between 3 x 1O~^M0 
and 6 x 10"''' M©. Our values for the hydrogen number 
density of between 3.5 x 10^^ cm~'^ and 9 x 10^^ cm""^, 
using a hydrogen atomic mass of 1.67x 10"^'* g, are equiv- 
alent to mass densities of ^ 10~^^gcm~^. Such values 
are consistent with the hydr ogen density derived for the 
Betelgeuse layer modeled bv IPerrin et al.l (j2007D . where 
they derived their hydrogen density by assuming a ra- 
tio between the II2O number density and the hydrogen 
numb er density of nn^o/nH ^ 10^^ (jJenning s & Sada 
[T99I and using their calculated values of nz/jO for the 
layer. Our values for the hydrogen density are there- 
fore also c onsistent with the ya lues for the water density 
derived bv IPerrin et al.l (|2007() , though they result from 
very different methods of calculation. 

4. THE ORIGIN OF THE POINT SOURCES IN THE MID-IR 

We now consider a physical interpretation of the point 
sources used to model the ISI interferometric data on 
Betelgeuse. Figure 4 shows images of the intensity dis- 



TABLE 4 

Number density and mass of the layer. 



Year 


nu (cm~^ X lO-"^-"^) 


Mlayer (Mq X lO"*) 


2006 


5±1 


6±1 


2007 


4.6±0.9 


3±1 


2008 


9±2 


4±2 


2009 


4.3±0.8 


3±1 


2010 


3.5±0.4 


5.0±0.9 



tribution of Betelgeuse during each observing epoch, cal- 
culated from the photosphere and layer models, with 
the point source locations indicated by filled white disks. 
The color brightness of the images is consistent with the 
same flux density scale, also depicted in Figure 4, and the 
areas of the disks representing the point sources indicate 
their flux densities calculated from their fltted fractions 
of the total intensity. All point sources are positive in 
intensity with respect to the star. 

Non-uniformities on the observed stellar disk of Betel- 
geuse at optical and near-IR wa velengths have been 
repor ted by several authors (e.g., iBuscher et al.l IT990l : 
Hau bois et al.lf2009| ). These asymmetries are generally 
modeled either as point sources, or as UDs or two- 
dimensional Gaussian proflles of fixed, unresolved size. 
The intensities of these 'spots' are generally found to be 
between 5% and 20% of the intensity of the stellar disk. 
The distribution and strengths of the spots are found to 
vary significantly across optical and near-IR band s dur- 
ing contemporaneous observations (jYoung et al.l [2000j . 
with the star sometimes appearing featureless in one 
band and highly non-uniform in another. The spots are 
also found to vary on timescales of less than 8 weeks 
(jWilson et al.lll997|). On e such spot was recently resolved 
in the ff-band bv [Haubo is ct al. (2009), who reported a 
diameter of O.Olf'for this spot. 

These resul t s are consistent with a prediction by 
iSchwarzschildl (|1975D for extremely large-scale convec- 
tion on red giants. Red supergiants are suggested to 
have only a few tens of convection cells on their sur- 
faces, as opposed to the thousands present on the surface 
of the Sun. Similar predictions also resu lt from three- 
dimensional numerical simu lations (e.g., iFrevtag et al.l 
I2OO2I : iChiavassa et al.l I2009D , and an initial attempt at 
matching simulations to in t erfero m etric observations was 
made by IChiavassa et al.l (|2010( ). ISchwarzschildl (|1975D 
predicted temperature fluctuations of up to ±1000 K 
with respect to the photospheric temperatures of red gi- 
ant stars. This is consistent with the intensities of spots 
observed in the visible and near-IR. 

The point sources fitted to our mid-IR data do not 
arise directly from large-scale convective features on the 
surface of Betelgeuse. If we assume a 0.01" diameter for 
the weakest of our point sources (see Table 2), which con- 
tributes 0.7% of the 4400 Jy flux density of Betelgeuse in 
2009, its temperature of ^660 K above the effective sur- 
face temperature, calculated using Equation 2, is consis- 
tent with a photospheric hotspot of ^1000 K behind the 
r = 1 layer. If we assign a similar diameter to any of 
the other point sources we observe, even greater temper- 
atures would result for associated photospheric hotspots. 
The mid-IR point sources are hence more intense than 
expected from surface temperature fluctuations caused 
by giant convection cells (Scliwarzschild . 1 975. ) . Further- 
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2006 - UD & PT ^^^H 2007 - UD 




-0.02 0.00 0.02 -0.02 0.00 0.02 

(arcsec) d^^ (arcsec) 



2008 - UD ^^^H 2009 - UD Sc 2 PTS 
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9p;^ (arcsec) 

Fig. 4. — Representations of the photosphere-and-layer model intensity distributions fitted to ISI observations of Betelgeuse for each 
observing epoch in the present work. The locations of the fitted point sources are shown as white disks: the areas of the disks scale linearly 
with their intensities. The visibility model fitted to each dataset is specified at the top of each imago: the 2006 data were fit with a UD 
and single point source, the 2007 and 2008 data were only fit with a UD, and the 2009 and 2010 data were fit with a UD and two point 
sources. The color bar shows the flux density per pixel for each color depicted, and each pixel is 1.2x1.2 milli-arcseconds in size. 
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more, every one of our point sources is located beyond 
the edge of the photospheric disk (see Figure 4). 

It is clear therefore that the observed asymmetry of 
Betelgeuse in the mid-IR is caused by a non-uniform in- 
tensity distribution in the layer. We do not have suffi- 
cient sampling of the UV plane in any epoch to identify 
whether the layer is in fact spherically symmetric, but 
with temperature fluctuations, or whether the layer is 
extended in particular directions. The point sources can 
also be explained in terms of opacity inhomogeneities in 
the layer. In the Rayleigh- Jeans approximation to the 
Planck law, the effective radiation temperature, Tb, of a 
column of material with optical depth r and temperature 
T is given by 

Ts=r(l-e-"). (10) 

For example, increasing the value of r from 1 to 3 in- 
creases Tb by a factor of 1.5. Non-uniform spatial, tem- 
perature and opacity distributions of the Betelgeuse layer 
could all explain the observed mid-IR non-uniformities. 

5. DISCUSSION 
5.1. The extended atmosphere of Betelgeuse 

Molecular layers directly surrounding the photospheres 
of red supergiant stars, and in particular Betelgeuse, have 
been modeled by a variety of authors under differing as- 
sumptions for the layer geometry and the dominant opac- 
ity source. A layer of water with a column density of 
10^°cm~^ at a temperature of ~1500K was first sug- 
gested by [TsuuI (f2000a) for Betel geuse, based on mod- 
eling of near-IR spectra. A combination of spectral and 
interferometric modeling by lOhna^ ()2004[ ). assuming a 
layer geometry identical to ours, suggested a layer tem- 
perature of 2050 K and an outer radius of 1.45i?*, while 
modeling the mid-IR opacity as being dominated by a 
continuum of water emission and absorption lines. Un- 
der th e assumption of a layer with a well-defined inner 
radius. pTsuiil (|2006[ ) modeled a shell beyond l.Si?* with a 
temperature of 2250 K, again dominated by water opac- 
ity. M ore advanced spectral modeling by Vcrhoclst et al.l 
([2006') revealed that water could not play a dominant 
role in determining the mid-IR continuum opacity, and 
instead suggested the presence of amorphous alumina in 
the molecula r layer. This model was consolidated by 
iPerrin et al.l ([2007) to suggest a shell between 1.32i?, 
and 1.42i?, at a temperature of 1520 K. 

The properties we derive for a layer surrounding Betel- 
geuse, summarised in Table 3 for each observing epoch, 
are largely consistent with these results. We find temper- 
atures that vary between 1900 K and 2800 K and outer 
radii that vary between 1.14i?, and l.SGR^,. Our model- 
ing te chnique most closely matched that of IPerrin et al.l 
(HOO?) in assuming thermal continuum emission from the 
layer within our observing band. This assumption is jus- 
tified by the featureless spectrum of Betelgeuse within, 
and around, our observing band and by the lack of any 
significant water line s within and surrounding our band 
(jWeiner et al.l[2003bD . The results are consistent with 
previous models of the Betelgeuse layer. 

A consistent picture of the composition of the layer has, 
however, not been achieved. While spectral modeling has 
uncovered the presence of numerous molecular species, 
including I I2O, C O, SiO and possibly CN, as well as alu- 
mma dust (frsuii|[2006- .Per rin et al.,,2007; .Kervella et al.l 



120091) ■ the opacity of the layer as a function of wave- 
length is not yet fully determined. Our results for the 
optical depth of the layer of ~1 at 1 1.15 iim wav e length 
are not in agreement w ith those of IPerrin et al.l ()2004l ) 
and IPerrin et al.l (|2007l ). who obtain optical depths of 
-2. 

The number density of atomic hydrogen, based on the 
11.15 /im opacity arising from electron- hydrogen atom 
collisions, is ~10"'^"'^ cm~^. This density, given the de- 
rived layer size, corresponds to H2O column densities 
of 10^^cm~^; this res ult is consistent with H2O column 
densities measured bv lJennings fc Sad a ( 1998) and those 
derived from spectral modeling (e.g., Tsuii 2006i) . We 
therefore suggest that electron-hydrogen atom collisions 
play a significant role in determining the mid-IR contin- 
uum opacity of the Betelgeuse layer, and that this mech- 
anism must also be accounted for in models for the layer 
composition. 

We have shown that the layer of material surrounding 
the photosphere of Betelgeuse has a non-uniform inten- 
sit y distribution. A sim ilar result was recently reported 
by lOhnaka et al.l (|2011[ ) based on interferometric mea- 
surements of the one-dimensional intensity and velocity 
profile of CO g as in the inner atmosphere of Betelgeuse. 
lOhnaka et al.l (|201lD also reported a change in the one- 
dimensional CO intensity and velocity profile between 
two epochs spaced by a year. We find that the inten- 
sity distribution of the layer varies between each the five 
observing epochs we present, which are also spaced by 
approximately one year. 

The non-uniformities that we observe in the layer re- 
flect, on much smaller spatial scales, the non-uniformities 
observed in various components of the circumstellar en- 
vironment of Betelgeuse. A gaseous component at sim- 
ilar temperatures to the Betelgeuse layer was show n to 
have an asymmetric intensity distribu tion by |l 7im et alj 
(|1998D . Our observations and those of iLim et al.i (| 19981 ) 
could be linked by the opacities in both wavelength 
bands being partly caused by electron-hydrogen atom 
collisions. The opacity due to this mechanism is pro- 
portional to implying that the observed stellar 
size decreases with shorter observing wavelengths. Ma- 
terial at chromospheric temperatu res is also observed 
to be asymmetrica l ly dis tributed (jHebden et al.l 119871 : 
iGilliland &: Duprei I1996D. and observation s of a cir- 
cumstellar envelope bv 'Ker vella "eTall ([20091) at near-IR 
wavelengths also reveal significant asymmetry. Asymme- 
try in the large-scale (~1") dust distribution surroundin g 
Betelgeuse was also observed by iBloemhof et al.l (|1985D . 

Molecular layer models have been applied to the few 
red supergiant stars that have been observed with suffi- 
cient spatial reso lution, such as u Cephei (Pcrri n et al.l 
[2OO5I iTsujil (20061 ) and a Herculis (jOhnaka 2 0041). Simi- 
lar layers are identified aro und AGB stars (jWeiner et al.l 
I2003at lOhnaka et aIll2005D . The properties of individ- 
ual red supergiant layers vary significantly from star to 
star, and a greater sample of measurements is required to 
understand what determines individual layer properties. 

5.2. Implications of the results reported for red 
supergiant mass-loss 

The processes by which red supergiant stars form and 
accelerate dust, while crucial to mass-loss scenarios, are 
little understood. Existing models of stellar atmospheres 
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Fig. 5. — Plot of the visual magnitude, and total 11.15 /xm flux 
density, of Betelgeuse over the period 2006—2010. The visual mag- 
nitude measurements, joined by a solid line, are 30-day averages 
of data obtained from the AAVSO online archive. The 11.15 ^m 
flux density measurements, also listed in Table 1, are shown here 
as crosses. The timespans of each ISI observing epoch are bounded 
by vertical dotted lines. 



(e.g.. iGustafsson et all [20081 iPled 12010( 1 cannot predict 
locations for dust condensation, and attempts to model 
dust acceleration through ra diation pressu re have also 
so far been unsuccessful fe.g.. lWoitkell2007l) . Our multi- 
epoch study of a layer directly surrounding Betelgeuse 
provides new insight into the processes involved in red 
supergiant mass- loss. 

We have shown that the basic properties of the layer, 
given in Table 3, exhibit striking variability from year 
to year. The layer is found to decrease in size and in- 
crease in optical depth over the years 2006 — 2008, and 
increase in size and decrease in optical depth between 
2009 — 2010. The layer temperature does not show such 
systematic variability, although the largest layer size (in 
2010) does correspond to the highest temperature. The 
variability in the ap parent mid-IR size of Betelgeuse re- 
ported bv lTownes et al. (2009.) can also be interpreted to 
indicate variations in the lay ers. Intriguin gly, a similar 
phenomenon was reported bv iPeas j ()1922l ). who noted a 
decrease in the size of Betelgeuse at visible wavelengths 
followed by an increase. 

Figure 5 shows measurements of the visual magnitude 
of Betelgeuse during between 2006 and 2010, as well as 
the ISI measurements of the total 11.15/im flux density 
(also listed in Table 1). The timespans of the ISI ob- 
serving epochs are also plotted. The visual magnitude 
measurements were obtained from the online databaseQ 
of the American Association of Variable Star Observers 
(AAVSO). The data, collected by a multitude of vol- 
unteer observers, were averaged over 30-day timespans; 
each timespan with less than 10 measurements was dis- 
carded. The visual variability of Betelgeuse over our ob- 
serving epochs was approximately 0.9 magnitudes. The 
visual intensity of Betelgeuse is not likely to represent 
that of the photosphere because the apparent stellar di- 
ameter at visible wa velengths is greater than that at 
near-IR wavelengths (jYoung et al.l l2000[ ). Also, using 

^ http:/ /www. aavso.org/data-download 



contemporaneous interferome tric observations of Betel- 
geuse at different wavelengths. I Young et al.l ()2000D found 
that hotspots evident in the visible were not present in 
the near-lR. Observations of Betelgeuse at visible wave- 
lengths could represent regions within the layer observed 
at mid-IR wavelengths. It is interesting that the mid- 
IR flux density of Betelgeuse, as well as the effective 
temperatures plotted in Figure 3, and the visual mag- 
nitude appear somewhat correlated, particularly in their 
increase during the 2009 and 2010 epochs. Previous au- 
thors (Kiss et al. 2006) have concluded that the visual 
variability of Betelgeuse is probably indicative of quasi- 
periodic pulsations, and is too great to be explained by 
variations in the hotspot distribution. 

The layer of material surrounding Betelgeuse is a can- 
didate location for the onset of dust formation. A num- 
ber of authors have identified the pres ence of alumina 
dust directly above the photosphere (iVerhoelst et al.l 
120061: IPerrin et al.l l2007t IVerhoelst et al.l 120091) . The 
derived physical parameters of the layer, such as the 
temperatures and hydrogen mass densitie s, are suitable 
for the formation of amorphous alumina (jWoitkd 120061 : 
iDirks et"al]l2008[ ). While some aspects of the modeled 
variability of the layer might be due to variations in the 
photosphere, measurements of changes in the apparent 
mid-IR appearance of Betelgeuse must be largely caused 
by variations in the layer. If alumina contributes sig- 
nificantly to the mid-IR continuum opacity, these vari- 
ations would represent variations in the dust content of 
the layer. Our estimates of changes in the mass of the 
layer, on the order of 10~^Mq (see Table 4), might indi- 
cate evolution in the layer. 

The point sources we observe in the layer of mate- 
rial surrounding Betelgeuse and their changes suggest 
that the layer dynamics are anisotropic. Asymmetries 
observed in the circumstellar environment of Betelgeuse 
are generally linked to the presence of giant convection 
cells on the stellar surface. Such convection cells are 
thought to shape the circumstellar environme nt of Betel- 
geuse by elevating cool photospheric material (iLim et aD 
[1998) and dri ving large- scale chromospheric motions 
([GiUiland fc Dupreelll996( ). Our results indicate that gi- 
ant convection cells have a significant role in shaping the 
Betelgeuse layer. This is because the mechanism that el- 
evates the layer material above the photosphere is clearly 
anisotropic. The large intensities of the observed point- 
source components imply the action of large-scale photo- 
spheric features, and their variability between observing 
epochs demonstrates the transient nature of the cause of 
the point sources. Changes in the distribution of con- 
vection cells on the photosphere of Betelgeuse could also 
infiuence the overall layer properties. While the nature 
of the non-uniform brightness distribution of the layer 
is yet to be determined, our observations represent evi- 
dence for the role of giant convection cells in shaping the 
inner atmosphere of Betelgeuse. 

5.3. Future work 

Observations of a sample of red supergiant stars with 
high spatial resolution and wide bandwidth, over multi- 
ple epochs, are necessary to form a general picture for 
the role of molecular layers in mass-loss processes. More 
detailed interferometric imaging of Betelgeuse at mid-IR 
wavelengths is also required to understand the true role 
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of giant convection cells in shaping the layer. Simulta- 
neous visible, near-IR, and mid-IR interferometry could 
also provide interesting information on correspondences 
between different regions of the layer. Interferometric 
studies with both high spatial resolution and hig h spec- 
tral resolution are planned for the ISI [Wishnow et al.l 
[2010), an d have been p erfor med at near-IR wave lengths 
bv .Ohnaka et al.l (|2009D and lOhnaka et al.l ()201lD . Such 
investigations are aimed at revealing the molecular com- 
positions, kinematic structures and true extents of red 
supergiant layers. More detailed modeling of layer 
opacities is also required, taking into account contribu- 
tions from molecular absorption and emission, dust and 
electron-hydrogen atom collisions. Upcoming optical and 
infrared interferometer systems, such as the Magdalena 
Ridge Observatory Interferometer and the MATISSE in- 
strument at the Very Large Telescope Interferometer, 
as well as high spatial resolution observations with the 
next generation of extremely large telescopes could fi- 
nally clarify our understanding of red supergiant mass 
loss. 

6. CONCLUSIONS 

The ISI has been used to conduct observations of Betel- 
geuse over 2006—2010 in the mid-IR continuum with high 
spatial resolution. By fitting visibility and closure phase 
measurements with image models, the apparent stellar 
size was measured at each observing epoch, and up to 
two point sources were sometimes detected at the edge 
of the stellar disk. The results were interpreted in terms 
of an absorbing layer surrounding an assumed 3641 K 
photosphere. We have shown that: 



1. Observations during every epoch can be modeled 
in terms of a layer surrounding the photosphere, 
with optical depths of '-^l at 11.15 /xm wavelength, 
outer sizes between 1.16i?» and 1.36-R*, and tem- 
peratures between 1900 K and 2800 K. 

2. Electron- hydrogen atom collisions contribute to the 
opacity of the layer. 

3. Assuming a non- varying photosphere, the layer ex- 
hibits significant variability. 

4. The layer has a non-uniform intensity distribution 
that also varies from year to year. 

Both the non-uniformity and the variability of the mod- 
eled layer suggest that giant convection cells play an im- 
portant role in shaping the inner atmosphere and the 
mass-loss dynamics of Betelgeuse. 
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